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Role of Spin-Orbit Coupling on the Spin Triplet Pairing in NaxCo02 ■ 2/H2O 
I: d-vector under Zero Magnetic Field 
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The d-vector in possibile spin triplet superconductor NaxCo02 • 3/H2O is microscopically 
investigated on the basis of the multi-orbital Hubbard model including the atomic spin-orbit 
coupling. As a result of the perturbation theory, we obtain the stable spin triplet supercon- 
ductivity where the p-wave and /-wave states can be stabilized. If we neglect the spin-orbit 
coupling, superconducting state has 6-fold (3-fold) degeneracy in the p-wave (/-wave) state. 
This degeneracy is lifted by the spin-orbit coupling. We determine the d-vector within the lin- 
earlized Dyson-Gorkov equation. It is shown that the d-vector is always along the plane when 
the pairing symmetry is p-wave, while it depends on the parameters in case of the /-wave state. 
The lifting of degeneracy is significant in the p-wave state while it is very small in the /-wave 
state. This is because the first order term with respect to the spin-orbit coupling is effective in 
the former case, while it is ineffective in the latter case. The consistency of these results with 
NMR and /iSR measurements are discussed. 

KEYWORDS: NaxCo02 ■ J/H2O; spin triplet superconductivity; multi-orbital model; d-vector 
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1. Introduction 

The unconventional superconductivity in strongly cor- 
related systems has been one of the central issues in the 
condensed matter physics. For example, heavy fermion 
superconductors, high-Tc cuprates,^^ organic supercon- 
ductors'^^ and Sr2Ru04^^ are cited. Recently, a new su- 
perconductor NaxCo02 -2/1120 was discovered by Takada 
et al,^^ and the possibility of unconventional supercon- 
ductivity has attracted huge interests. 

Immediately after the discovery of superconductivity 
in NaxCo02 • 2/H2O, the symmetry of superconductiv- 
ity has been studied by many experimental measure- 
ments. While some controversial results exist, 
most of them have indicated the non-s-wave super- 
conductivity. For example, the absence of coherence 
peak in NMR l/TiT^'^) and the power low behaviors 
in and specific heat^"^"^^-* are evidences for 

the anisotropic pairing. Recently, a magnetic phase has 
been discovered in the neighborhood of superconducting 
phase. This observation clearly indicates the impor- 
tance of electron correlation which generally leads to the 
non-s-wave superconductivity. 

These experimental indications have accelerated the- 
oretical studies on the superconductivity in NaxCo02 • 
J/H2O. In the first stage, the effect of frustration stimu- 
lated many interests since NaxCo02 • 2/H20 has layered 
structure constructed by the triangular lattice. The RVB 
theory has been applied to the triangular lattice^^"^'') 
and concluded the spin singlet d-wave superconductivity. 
Some authors have pointed out the frustration of charge 
ordering at the electron filling n = 4/3,^^-' and the /-wave 
superconductivity due to the charge fluctuation has been 
discussed. ^^'^^^ Recently, the RVB theory has been ap- 
plied to the multi-orbital model and concluded the spin 
triplet superconductivity.^^^ 

Another theoretical approach is based on the pertur- 
bation expansion from the weak coupling region which 



includes the perturbation theory,^^^ random phase ap- 
proximation, FLEX approximation^''^ and perturba- 
tive renormalization group. '^^^ Some authors have taken 
account of the Fermi surface of NaxCo02 • J/H2O partly, 
and concluded the /-wave superconductivity. 
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superconductivity,^^^ i-wave superconductivity"^^' and 
nearly degeneracy between d- and /-wave superconduc- 
tivities. 

Although these theories except for Refs. 27 and 29 have 
assumed single-orbital models, the conduction band in 
NaxCo02 • yH20 has orbital degeneracy, as pointed out 
by Koshibae et aZ."^^' The conduction band mainly con- 
sists of three t2g-orbitals in Co ions which hybridize with 
02p-orbitals. Therefore, it is highly desired that the pair- 
ing symmetry in this material is examined on the basis 
of the multi-orbital model. 

For this purpose, we have constructed a three-orbital 
Hubbard model which appropriately reproduces the elec- 
tronic structure obtained in the LDA calculations, 
and we have analyzed it on the basis of the FLEX 
approximation^^) and perturbation theory."^*' Then, we 
have obtained some notable results which are summa- 
rized as follows, (i) The spin triplet superconductivity 
is stable unless the Hund's rule coupling is very small, 
(ii) The p-wave state and /-wave state are nearly de- 
generate owing to the orbital degree of freedom, (iii) 
There is a nearly ferromagnetic spin correlation along 
the plane which stabilizes the spin triplet pairing. This 
is consistent with recent neutron scattering measure- 
ments which have reported the anti-ferromagnetic or- 
der with stacking ferromagnetic plane. "^^^ (iv) The ver- 
tex correction whose importance has been pointed out 
for Sr2Ru04'**'''*^' is not important in NaxCo02 • yB.20. 
(v) The two of three orbitals, namely Cg-doublet are es- 
sential for the superconductivity. The orbital-dependent- 
superconductivity proposed for Sr2Ru04^^' is partly jus- 
tified in NaxCo02 • 2/H2O. (vi) However, the orbital de- 
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generacy in Cg-doublet is particularly important because 
the single-orbital approximation artificially suppresses 
the p-wave superconductivity. This point is in sharp con- 
trast to Sr2Ru04 where the single-orbital approximation 
is valid. '^"^ In this sense, NaxCo02 • yH20 is a typical 
multi-orbital superconductor in d-electron system. 

One of the interesting subjects in the multi-orbital su- 
perconductor is the role of spin-orbit coupling. This issue 
is particularly important in the spin triplet superconduc- 
tivity which has an internal degree of freedom described 
by the d- vector. ^"'^ If we neglect the spin-orbit coupling 
as in the previous studies, the 6-fold (3-fold) degen- 
eracy remains in the p-wave (/-wave) state at T = T^. 
Therefore, we have to take into account the spin-orbit 
coupling to determine the pairing state. The goal of this 
paper is to clarify the role of spin-orbit coupling and mi- 
croscopically determine the d-vector in NaxCo02 • t/H20. 

The rf-vector is particularly important to discuss the 
Knight shift measurement which has been a powerful 
method to determine the pairing symmetry.^^^ This is 
because the magnetic susceptibility in spin triplet su- 
perconductor significantly depends on the direction of 
d-vector. Although many experiments have been per- 
formed to determine the pairing symmetry in NaxCo02 • 
j/H20,6-i*^) there is any conclusive evidence neither for 
the spin triplet pairing nor for the spin singlet pairing. 
We think that this is partly due to the lack of knowl- 
edges on the pairing state expected in the spin triplet 
superconductivity. The results in this paper will provide 
a clear subject for a comparison between the theory and 
experiment. 

The role of spin-orbit coupling on the spin triplet su- 
perconductivity has been a longstanding problem since 
the discovery of heavy fermion superconductors.'''*^^) 
However, the microscopic study has not been performed 
owing to the complicated electronic structure of heavy 
fermion systems. The discussion about the pairing sym- 
metry in UPts still continues'*''"^^) partly because there 
is no microscopic research on the anisotropy of rf- vector. 
Recently, we have developed a microscopic theory on the 
d- vector and applied to Sr2Ru04.^''' The present study 
on NaxCo02 • yll20 provides a contrasting example and 
we expect that these studies on the d-electron systems 
will lead to a systematic understanding including the /- 
electron systems. 

This paper is organized as follows. In §2, we introduce 
the three-orbital Hubbard model including the atomic 
spin-orbit coupling and derive the two-orbital Hubbard 
model like in the previous study. '^^^ The pairing symme- 
try allowed in this model is classified in §3. The linearized 
Dyson-Gorkov equation in the multi-orbital model in- 
cluding the spin-orbit coupling is developed in §4.1. In 
§4.2, the pairing state is determined on the basis of the 
second order perturbation theory. We show that the role 
of spin-orbit coupling is quite different between the p- 
wave superconductivity and /-wave one. This difference 
is illuminated by showing the splitting of in §4.3. We 
show that the rf-vector in the p-wave state is strongly 
fixed against the magnetic field, while that in the /-wave 
state is fixed weakly. In §4.4, we discuss the cross-over 
from the weak spin-orbit coupling region \ ^ W rele- 
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vant for NaxCo02 • y^20 to the strong spin-orbit cou- 
pling region A » VF. Although the latter is unrealistic 
for NaxCo02 • yH20, this analysis will be useful for a 
unified understanding including the heavy fermion su- 
perconductors. In §5, we summarize the comparisons be- 
tween NaxCo02 • yH20 and Sr2Ru04. Some discussions 
are given in the last section §6. 

2. Spin-Orbit Coupling in Multi-Orbital Hub- 
bard Model 

First of all, we introduce a three-orbital Hubbard 
model including the spin-orbit coupling from which a 
two-orbital model are derived later. We consider a two- 
dimensional model which represents Co ions on the trian- 
gular lattice. We have constructed a tight-binding model 
for Co t2g-orbitals which reproduces the results of LDA 
calculations.^^) By adding the spin-orbit coupling term 
and Coulomb interaction term, the three-orbital Hub- 
bard model is obtained as, 

Hs = Ho + Hj^s+Hi, (1) 

i,j,s a,b 

Hj^s = 2XY,LiSi, (3) 

i a i a>b 

i a>b 

Here, the indices i and j denote the sites in the real space 
and indices a and b denote the orbitals. We assign the 
dxy-, dyz- and rfxz-orbitals to a = 1, a = 2 and a = 3, 
respectively. 

The first term Hq is a tight-binding Hamiltonian where 
ta,bs,j are determined according to the symmetry of or- 
bitals and lattice. The dispersion relation expected in the 
LDA calculation is reproduced by assuming nine hopping 
parameters from ti to tg as well as the crystal field split- 
ting Cc which arises from the distortion of octahedron. 
For instance, we assume ti,i,i,i±a = ti, j±b ~ *2, 
i2,3,i,i±a = ^35 *l,l,i,,±(a-b) ~ *i.i,i.»±2a = Hi 

i2,3,i,i±2a = U, i2,3,i,i±(a-b) = *i,3,i,i±(a-b) = 
and 2 i j±(a-b) ~ ^® denote the basis of triangular 
lattice as a=(\/3/2, — 1/2) and b=(0, 1) and we choose 
the lattice constant as a unit length. The other hopping 
matrix elements ta,b,i,j are obtained by the symmetry op- 
eration. Since ts is largest among tn, we choose the unit 
of energy as ts = 1 throughout this paper. We describe 
Hq in the matrix representation as, 

k,s 

where cl = (cl ,c\ ,c\ ) is a vector representa- 

k,s ^ fe,l,s' fe,2,s' k,3,s' ^ 

tion of Fourier transformed creation operators with spin 
s. The matrix element of H{k) has been given in Refs. 37 



Full Paper 



J. Phys. Soc. Jpn. 



Full Paper 



Youichi Yanase, M. Motizuki and Masao Ogata 



3 



and 38. 

In order to clarify the nature of superconductivity in 
this model, it is useful to introduce a non-degenerate 
aig-orbital and doubly-degenerate Cg-orbitals. They are 
defined as, 

>=^(|xy>+|xz>+|yz>), (6) 



|eg,i >= -^(1^2 > -\^y >)' 



(7) 



|eg,2>=-^(2|yz>-|xz>-|xy>). (8) 

The orbital-dependent-superconductivity in NaxCo02 • 
yB.20 occurs in this basis as shown in Refs. 37 and 38. 
We choose a basis of Cg-orbitals different from Ref. 38 
in order to simplify the following notations. By choosing 
the basis wave function as eqs. (6-8), the tight-binding 
Hamiltonian is transformed to be, 



k 



and dl 
and 



(9) 



are 
,2 >- 



through the unitary transformation. Here, 
(d\ ,d\ ,d\ ) where d] , d\ 

creation operators of |aig >, jcg, 1 > 
orbitals, respectively. 

Wc choose the c-axis as a quantization axis of spin 
for a convenience of following discussion. By choosing 
the basis functions as eqs. (6-8), we obtain a simplified 
expression for the atomic spin-orbit coupling term Hi,s 
as, 
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.(10) 



Interestingly, the matrix element of eq. (10) is the same 
as that in Sr2Ru04.®°'^^^ Therefore, we can discuss the 
role of spin-orbit coupling in NaxCo02 • yH2 in anal- 
ogy with Sr2Ru04. Note that the basis wave function 
in Sr2Ru04 is d^y-, dyz- and dzx-orbitals, while that is 
Wig >7 kg, 1 > and |eg, 2 >-orbitals in NaxCo02 • JI/H2O. 
This difference arises from the position of apex oxy- 
gens. While an apex of RuOe octahedron is along the 
c-axis, all of apex oxygens in CoOg octahedron are tilted 
from the c-axis. When we consider the matrix element in 
eq. (10), the dxy-orbital in Sr2Ru04 corresponds to the 
aig-orbital in NaxCo02 • J/H2O, while the dyz- and d^x- 
orbitals correspond to the Cg-doublet. It is expected that 
the dxy-orbital is responsible for the superconductivity 
in Sr2Ru04,^'^'^'^'' while the superconductivity is mainly 
caused by the Cg-doublet in NaxCo02 -2/1120.3^' Thus, 
these two materials provide contrasting examples of spin 
triplet superconductors. 

The coupling constant 2A in f^LS has been estimated as 



57mev.^^^ This value corresponds to A = 0.17 in our unit 
if we choose the band width W ^ dt^ = 1.5eV according 
to the LDA calculations. Since this estimation has some 
ambiguities, we investigate the range from A = 0.05 to 
A = 0.25. In this range, the efi'ect of spin-orbit coupling 
on the band structure is very small. There is a hole pocket 
enclosing the F-point and six hole pockets near the K- 
points, which are consistent with LDA calculations. ^''"^^^ 
The former (oig-Fermi surface) mainly consists of the 
aig-orbital and the latter (cg-Fermi surface) mainly con- 
sists of the Bg-orbitals. 

From the analysis of three-orbital Hubbard model with 
A = O,^^'^^-* we have found that the superconductivity is 
mainly induced by the Cg-orbitals and the Oig-orbital can 
be simply ignored to discuss qualitative features of super- 
conductivity. Therefore, we derive a two-orbital Hubbard 
model by simply dropping the aig-orbital as. 



H2 = H, 



(2) 



r(2) 



fe,s 



(11) 

(12) 



i a=l i a>b 



i a>b 



i a^b 



(13) 



Here, we have introduced a two component vector = 
(dl , dl ) and 2x2 matrix 



Here, 



e[,{k) 
(fc) -|- isA 

is obtained from 



5'l2(fc) - 



(14) 



(9) 



as e 



,(fc) 



eq. 

H'{k)i+ij+i where H'{k)ij is the matrix element of 
H'{k). It should be noticed that the off-diagonal ma- 
trix elements connecting up and down spins in eq. (10) 
vanish. This means that the operators Lx and Ly have 
no matrix clement in the Hilbert space expanded by Cg- 
orbitals. Therefore, only the spin-orbit coupling along the 
>axis is effective in the two-orbital model. This fact re- 
markably simplifies the calculation and provides a clear 
understanding for the results in §4. 

The second term in eq. (11) describes the on-site 
Coulomb interactions including the intra-orbital repul- 
sion U, inter-orbital repulsion U', Hund's rule coupling 
Jh and pair hopping term J. We impose the relations 
U = U' + Jb_ + J and Jn = J throughout this paper. 
The Coulomb interaction term is invariant for the uni- 
tary transformation with respect to the orbital owing to 
these relations. 

In this paper, we investigate the d- vector in NaxCo02- 
J/H2O on the basis of the two-orbital Hubbard model in 
eq. (11). According to Refs. 37 and 38, we choose the 
hopping parameters as, 



4 J. Phys. Soc. Jpn. 



Full Paper 



Youichi Yanase, M. Motizuki and Masao Ogata 



a(0.1, 0.2, 0.3, -0.2, -0.05, 0.2, 0.2, -0.25). (15) 

where a = 0.6 ~ 1 is consistent with LDA calculations. 
We choose a small value of a = 0.6 in this paper. This is 
because the electric DOS in Cg-Fermi surface decreases 
by neglecting the oig-orbital and this artificial decrease 
is compensated by decreasing a. We define the number of 
holes in the two-orbital model as rie which corresponds to 
the area of eg-Fcrmi siirfacc. According to the chemical 
analysis, the total number of holes is n = 0.4 0.5 
which is different from n = 0.67 in the parent mate- 
rial Nao.33Co02. Since < no < n, wo investigate the 
range = 0.05 ~ 0.35. We have found that rie is an im- 
portant parameter for the superconductivity rather than 
38) rpjjg superconducting instability is significantly sup- 
pressed when rie = 0. For finite n^, the superconduct- 
ing Tc slightly increases with the increase of crystal field 
splitting Cc which leads to the increase of rie.^*'^^-' The 
recent NMR measurement has confirmed an important 
role of e^}^^ We show a typical Fermi surface of two- 
orbital model in Fig. 1. 



M 



)^0 OK 



Fig. 1. Fermi surface in the two-orbital Hubbard model (solid 
lines). We have shown the Fermi surface of aig-orbital which is 
determined as H'(k)ii — fi = (thin dashed line). The thick 

(laslicd line is the first Brillouin zone of triangular lattice. The 
parameters are chosen to be A = 0.17, a = 0.6, rie = 0.21 and 
n = 0.5. 



are two wave functions in both p-wave symmetry and /- 
wave symmetry. The Px-wave and Py-wave are degenerate 
owing to the symmetry of triangular lattice. On the other 
hand, the /i-wave and /2-wave are not degenerate. 



El 


Px-wave 


sin cos ^ky 


Py-wave 


^ — ^ — ^7: 

sin ky + sin ^ ky cos k^ 


Bi 


/i-wavc 


sin ^/!;y(cos ^ky — cos -^kx) 


B2 


/2-wave 


sin ^fex(cos |/cy — cos ^kx) 



Table 1. Orbital part of odd-parity pairing function in the tri- 
angular lattice. The first column shows the irreducible repre- 
sentations of Dg-group. The second column shows the notations 
adopted in this paper, which are the counterparts of the isotropic 
system. The third column shows the typical wave function. Note 
that the wave function obtained in the Dyson-Gorkov equation 
is different from the third column to some extent. 



The spin part of pairing function is described by the 
d- vector as,**'^'^'') 



AxT(fc) Au(fc) j -'^W'^^y 



(16) 



-d^{k) + idy{k) d^{k) 

d^{k) d^{k)+idy{k) 

Since the Cooper pair has spin 5=1, there is a 3- fold 
degeneracy in the SU{2) symmetric system. Therefore, 
if the spin-orbit coupling is neglected, the degeneracy in 
the p-wave state is 3 x 2 = 6-fold due to the spin part 
and orbital part, while that in the /-wave state is 3-fold. 





Px£ + PyV, PyX - Pxy 


fxy — 


Px£ - PyV, PyX + p^y 


Pz 


(Px ±iPy)z 




fix - af2y, af2X + fiy 




hi 



Table II. Classification of the pairing state including the spin- 
orbit coupling. The first column shows the notations adopted in 
this paper. 



It should be noted that ARPES measurements^^' for 
non-superconducting Naa;Co02 observed the oig-Fermi 
surface, but the Cg-Fermi surface has not been found. 
This observation implies rie = which contradicts with 
our basic assumption, namely a finite value of rig. Our 
analysis on the three-orbital Hubbard model has shown 
that the superconductivity is hardly stabilized when 
rie = 0.^^) This result implies that the Cg-Fermi sur- 
face exists in the superconducting materials as expected 
in LDA calculations. The increase of Ue is actually ex- 
pected in superconducting materials because the H2O 
molecules increase the crystal field splitting Cc, as shown 
by the NMR measurement.^®^ 

3. Classification of the Pairing State 

Before the microscopic analysis, wc classify the sym- 
metry of pairing state in the two-orbital Hubbard model 
including the spin-orbit coupling. First, we show the or- 
bital part of spin triplet pairing function in Table I. There 



Table II shows the classification of pairing state for 
finite spin-orbit coupling. According to the results for 
A = 0,^^^ the /2-wave state is not stabilized in the multi- 
orbital Hubbard model. Therefore, we have not shown 
the pairing state which is reduced to the /2-wave state 
in the limit A — > 0. Since the spin-orbit coupling A < 0.25 
is much smaller than the band width W ^ \a\ in Fxy- 
state is much smaller than 1. 

The Fxy-state and all of the p-wave states in Table II 
are two-dimensional representations. Although the two 
pairing states in Pxy-i- are not degenerate in general, "^^^ 
there remains an additional degeneracy in the two-orbital 
Hubbard model. This is because the matrix element of 
and Ly vanishes in this model. The symmetry of Hubbard 
Hamiltonian with A = is G SU{2) T ® C/(l), where 
G, SU{2), T and U{1) show the symmetries of point 
group, spin rotation, time-reversal and gauge transfor- 
mation, respectively. When the spin-orbit coupling ex- 
ists, the symmetry is reduced to G U{1) (g) T U{1). 
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Thus, the SU{2) symmetry is violated by the spin-orbit 
coupling. However, the U{1) symmetry for the spin ro- 
tation in the plane remains because = Ly = 0. This 
additional U{1) symmetry is the origin of degeneracy in 
Pxy4--state. In the three-orbital Hubbard model (eq. (4)), 
the symmetry is reduced to G^T® [/(I), and the degen- 
eracy in the Pxy+-state is lifted. However, it is expected 
that the lifting of degeneracy is small because the order 
parameter in oig-orbital is very small. 

4. Determination of the Pairing State 

4-1 Linearized Dyson-Gorkov equation without SU{2) 
symmetry 

In this section, we determine the pairing state in the 
two-orbital Hubbard model by solving the linearized 
Dyson-Gorkov equation within the second order pertur- 
bation theory. The derivation of Dyson-Gorkov equation 
has been explained in literatures. The application to 
the multi-orbital model including the spin-orbit coupling 
has been given in the study of Sr2Ru04.^"^ The applica- 
tion to the two-orbital model (eq. (11)) is simpler than 
that to Sr2Ru04 because the transverse component of 
spin-orbit coupling vanishes. 

In order to make following discussions clear, we intro- 
duce a unitary matrix U{k, s) = {uij{k, s)) which diago- 
nalizes H^^\k,s), namely. 



U^k,s)H^'^\k,s)U{k,s) 



Ei{k) 






E2{k) 



(17) 



Here, Ei(k) < E2{k) and Ei{k) do not depend on the 
spin owing to the time-reversal symmetry and inver- 
sion symmetry. The Green function G{k,s) = (it<;„l — 
H^'^\k,s))~^ is obtained as, 

2 

Gij{k,s) = ^Uiaik,s)u*„{k,s)Ga{k), (18) 

a=l 

where Ga{k) 



(a) 

a,si^ 



(b) 



(c) 



b,s2> 



■*-C,s3 



7-*- 



■»-d,s4 

(d) 







> •• '•') 



(e) 



->Tr 



7^ 



Fig. 2. (a) Diagrammatic representation of the irreducible vertex 
in the particle-particle channel, (b-e) The second order terms 
with respect to the Coulomb interactions (dashed lines). The 
solid line denotes the Green function having the indices of spin 
and orbital. 



equation is written in terms of an effective interaction 
within the i?2(fe)-band. 



{k) = - Vik, k', si, s2, s3, s4) 

fe',s3,s4 

x|G2(fc')|'A,3.4(fc'), (19) 



with 



V{k, k', si, s2, s3, s4) = J2 <2(fc, sl)<2(-fc, s2) 



abed 



xVabcdik, k', si, s2, s3, s4)uc2(fe', s3)Md2(-fe', s4). (20) 

Here, Vabcd{k, k', si, s2, s3, s4) is the irreducible vertex in 
the particle-particle channel having the indices of orbital 
and spin (sec Fig. 2(a)). Because of = Ly = 0, the z- 
componcnt of spin is conserved and V{k, k' , si, s2, s3, s4) 
is finite only if si -|- s2 = s3 + s4. Therefore, the lin- 
earized Dyson-Gorkov equation is reduced to the eigen- 
value equations for An, A^j^ and Aj^j^, respectively. 

A|TAn(fc) = -^T/n(fc,fc')|G2(fc')pATT(fc'), (21) 

k' 

Al^An(fc) = -5]Fn(fc,fc')|G2(fc')l'An(fc'), (22) 

k' 

Ai^Au(fc) = -l]ni(^,fc')|G2(fc')l'Au(fc'), (23) 

k' 

where V^^{k,k') - ^(fc. A:', T, T, T, T), Vi^{k,k') = 
y(fc,fc',i,i,i,i) and Vn{k,k') = V {k,k', ],[,],[) + 
V{k, k', t, i, i, t)- It is notable that the maximum eigen- 
values of eqs. (21) and (23) are equivalent. The wave 
functions are related to be An^k) = e"'^A|^(/c) where 
the phase 4> is arbitrary. 

If we neglect the spin-orbit coupling, the maximum 
eigenvalue does not depend on the direction of c?-vector 
since V^^{k,k') — Vii{k,k') — V'^i{k,k'). However, we 
find A^ = X^^ ^ Xl^ when A is finite. The super- 
conducting transition is determined by the criterion, 
max(A|T,Al^) = 1. If Ap > A^^, the 
Fxy-state is stabilized at T = Tc, while P^- or F^-state is 
stabilized when A]^ < XlK 

In this paper, we estimate irreducible vertex 
Vabcd{k, k', sl, s2, s3, s4) up to the second order with re- 

(2) 

spcct to the Coulomb interaction term i/j . The dia- 
grammatic representation is shown in Figs. 2 (b-e) which 
are calculated from the possible combination of Coulomb 
interactions and Green functions. We numerically solve 
the eigenvalue equations for A|^ and A]-*- and determine 
the pairing state with highest T^. In the following, we di- 
vide the first Brillouin zone into 128 x 128 lattice points 
and we take 512 Matsubara frequencies for T > 0.005 
and 1024 Matsubara frequencies for 0.002 < T < 0.005. 
We have confirmed that this size of calculation is suffi- 
cient for the following results. 



As shown in the LDA calculations, the energy band 
described by E2{k) crosses the Fermi level and Ei{k) 
is far below the Fermi level. In this case, the supercon- 
ducting transition is induced by the Cooper pairing in 
the E2 (fe)-band. Therefore, the linearized Dyson-Gorkov 



4-.2 d-vector 

First, we show the phase diagram of two-orbital Hub- 
bard model without spin-orbit coupling which will be a 
basis of following discussions. Figure 3 shows the stable 
pairing symmetry in the phase diagram of and Jh- In 
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the figure, the intra-orbital repulsion U is determined so 
as to obtain Tc = 0.01 which corresponds to 20K in our 
unit. This temperature is higher than the observed tran- 
sition temperature Tc ^ 5K. However, the stable pairing 
symmetry is almost independent of the temperature as 
shown in Ref. 38. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 

n 



Fig. 3. Phase diagram without spin-orbit coupUng in the nc-Jn 
plane. The soUd Une is the phase boundary obtained by the in- 
terpolation. 



Among many parameters in the multi-orbital Hubbard 
model, the relevant parameter for the pairing symme- 
try is the Hund's rule coupling Jh and the number of 
holes in eg-Fermi surface nc.'^*' Figure 3 is quite simi- 
lar to the phase diagram of three-orbital model shown in 
Ref. 38. When the Hund's rule coupling is large and/or 
Tie is small, the p-wave superconductivity is stabilized. 
The /-wave superconductivity is stabilized for small Jh 
and/or large rio. The only qualitative difference between 
the two-orbital and three-orbital models is that the /- 
wave state stabilized for Uc < 0.1 and J-a/U > 0.2 in 
Ref. 38 vanishes in Fig. 3. 
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Fig. 4. Phase diagram including the spin-orbit coupling A = 0.17. 
The notations of Pxy+ , ^xy and have been given in Table II. 
The solid line is the phase boundary obtained by the interpola- 
tion. 



We show the results of d-vector in Fig. 4, where the 
Tc is fixed to be Tc = 0.003 ~ 6K consistent with exper- 
imental value. Hereafter, the parameter U is determined 
so as to obtain this Tc. For example, we obtain U =^ 6.57, 
U = 5.47 and U = 4.34 in A, B and C in Fig. 4, respec- 
tively. It is clearly shown that the spin-orbit coupling 
stabilizes the Pxy+-statc when the pairing symmetry is 
p-wave. Then, the d-vector is parallel to the plane. On 
the other hand, the d-vector in the /-wave symmetry de- 
pends on the parameters, although the Txy-state seems 
to be most stable. 

In order to understand these results, it is useful to 
analyze the role of spin-orbit coupling in a perturbative 
way. Although we have included the spin-orbit coupling 
term into the unperturbed Hamiltonian, the perturbative 
treatment is valid since A is much smaller than the band 
width. As a result of perturbation expansion in A, the 
effective interaction Vss'ik,k') is written as. 



OG 



Vss' {k, k') = [k, k') + k') 

n=l 



(24) 



We have actually applied this perturbation theory to 
Sr2Ru04 and shown that the first order term vanishes 
when the 7-band is mainly superconducting. ^''^ This is 
because the hybridization term disappears between dxy- 
orbital and dy^- dzx-orbitals. However, the first order 
term exists in case of NaxCo02 • yH20 since the large 
hybridization term exists in the Cg-doublet. 

It is easy to find that the first order term always ap- 
pears with the combination to the off-diagonal Green 
function or off-diagonal matrix element of U{k,s). Tak- 
ing into account the symmetry of e'i2i^) ^'i- (14), 
namely e[2{k^,ky) -ei2(-fcx, fcy) = -e^^ik^,, -ky), we 
obtain the following relations. 



vll\k,k')^vll^*{k,k'), 

V^^\k:xT ky, ky,, ky) = —Vj^\—kx,ky,—k^,ky) 



(25) 



-Vj^^k-^, fcy,/c^, fcy), (26) 



Vll\k,k') = 0. 



(27) 



According to eq. (26), the kernel of linearized Dyson- 
Gorkov equation works on the p^-wave state to produce 
the py-wave state, namely 

^Ay(fc)l/«(fc,fc')|G2(fc')l'Ax(fc') ^0, (28) 



and 



^Ax(fc)KV)(fc,fc')|G2(fc')l'A, 



k,k' 



ik') 



j2Ay{k)Vll\k,k')\G2{k')\'Ay{k')^0, 

k,k' 



(29) 



where Ax(fc) and Ay(/c) are functions having the sym- 
metry of px-wave and Py-wave, respectively. Thus, the 
first order term in A works like a first order perturba- 
tion on the degenerate states, namely px- and py-states. 
Therefore, the eigenvalue Aj^ for one of the T'xy+- and 
Pyiy — -states increases in the first order of A, while 
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for the other state decreases. On the other hand, the 
first order term vanishes for the d-vector along the c- 
axis according to eq. (27). Therefore, the rf- vector along 
the plane is always stabilized within the first order the- 
ory in A independent of the microscopic details. Which 
state is stabilized between Pxy+ and Pxy- is determined 
by the microscopic calculation for Vss'{k,k'). Figure 4 
shows that the Pxy+-state is stable within the second 
order perturbation theory. 

The situation is quite different for the /-wave super- 
conductivity. Owing to eq. (26), the first order term is 
not eff'ective in the /i-wave state, namely 



Y,^ii{k)vff{k,k')\G^ik')\'An{k') = 0, 



(30) 



k,k' 



where Afi(fc) is the function having the; /i-wave sym- 
metry. Although the first order term couples the /i-wave 
state to the /2-wave state, this is a higher order effect be- 
cause; thc;sc states are non-degenerate. Indeed, the eigen- 
value of linearized Dyson-Gorkov equation for the /2- 
wave state is much smaller than that for the /i-wave 
state. Then, the role of first order term, V-I^\k,k') is 
negligibly small. Therefore, the direction of d-vector is 
determined by the higher order terms beyond the second 
order, namely v}"\k, k') (n > 2). The role of higher or- 
der terms depends on the microscopic details as we have 
shown in the case of Sr2Ru04.^°^ In NaxCoOz • yHaO, 
Fxy-state (i^z-state) is stabilized for small (large) value 
of rie and/or large (small) value of Jh as shown in Fig. 4. 



(a) P-wave 

first order p higher order 

xy- P 

■ — ^ "y- 

6-fold y' ^ P, ^ 

X Pxy^ Pxy. 

(b) F-wave 

3-fold 3-fold ...^ ^ 



Fig. 5. Schematic figure of level scheme for (a) p-wave supercon- 
ductivity and (b) /-wave superconductivity. The degeneracy at 
A = is shown on the left. The level splitting due to the first 

order perturbation in A is shown in the center. The energy lev- 
els affected by the higher order perturbation are shown on the 
right. We assume in (b) that the Fxy-state is stable rather than 
the Fz-state. 



We summarize the above discussions in Fig. 5. In case 
of the p-wave superconductivity, the first order pertur- 
bation in A stabilizes the Pxy-i--state and destabilizes the 



xy- 



-state. The P^-state is not affected in this order. We 
find that higher order perturbation terms favor the Pxy-i-- 
state furthermore, as shown in Fig. 5(a). In case of the 
/-wave superconductivity, three-fold degeneracy is not 



lifted by the first order perturbation. The higher order 
perturbation stabilizes the Pxy-state or P^-state depend- 
ing on the parameters such as Jh and rie- 

4-3 Splitting of Tc 

The quite different role of spin-orbit coupling between 
the p- and /-wave symmetries discussed above is illumi- 
nated by showing the splitting of Tc- In eqs. (21-23), the 
linearized Dyson-Gorkov equation is defined for each di- 
rection of rf-vector. We define the superconducting Tc for 
d \\xy and that for d \\zhy the criterion X^Tc^^'') = 1 
and Xl^{Tc^^) = 1, respectively. The former corresponds 
to the Tr of P, 



xy4-' 



or Pxy-state while the latter is the Tc 



of Pj;- or P^-state. 

Figure 6 shows the splitting of Tc, namely AT^/T^ 



(xy) 



■Ti^Vl^c 



The parameters Jh and ric are chosen 



to be A, B, C in Fig. 4, where the Pxy+-state, Pxy-state 
and Pz-state are stable, respectively. It is clearly shown 
that ATc is much smaller in the /-wave symmetry than 
in the p-wave symmetry. This is because the first order 
term in A is ineffective. Indeed, while ATc/Tc increases 
linearly in the p-wave symmetry, the increase in the /- 
wave symmetry is nearly proportional to the square of 
A. 



0.25 
0.2 



A, p-wave 

B, f-wave 

C, f-wave 




0.25 



T^^'>\/Tc. The pa- 



Fig. 6. A-dependence of ATc/Tc = |Tc 

rameters are chosen to be A, B and C in Fig. 4. The Pxy+-state 
Fxy-state and F^-state are stable, respectively. 



The quite different magnitude of ATc/Tc shown in 
Fig. 6 will be reflected in the phase diagram under the 
magnetic field. In order to gain the Zeeman coupling en- 
ergy, the d- vector can be rotated by the applied magnetic 
field. Fig. 6 indicates that the d-vector is strongly fixed 
against the magnetic field in case of the p-wave symme- 
try, however the d- vector in the /-wave symmetry rotates 
in a wc^ak magnetic field. We have actually determined 
the phase diagram under the magnetic field on the basis 
of the weak coupling theory and found that the rotation 
docs not occur in the p-wave symmetry up to the half of 
Pauli paramagnetic limit. '"'^^ The results for the multiple 
phase diagram under the magnetic field and the charac- 
teristics of each phase will be shown in the subsequent 
publication.^®) 



8 J. Phys. Soc. Jpn. 



Full Paper 



Youichi Yanase, M. Motizuki and Masao Ogata 



The resuhs of Figs. 4 and 6 are important for an inter- 
pretation of NMR and /xSR results because the magnetic 
susceptibility has an anisotropy arising from the direc- 
tion of d- vector. "^^'"'^^ In general, the magnetic suscepti- 
bility in the unitary state is obtained as, 

Xtiv/ Xn — 

< V - {di,{k)d.{k)/\d{k)\^){l - Y{k,T)) >F, (31) 

where Of is an average on the Fermi surface, < A >f= 
p(0)~^/ Ad{E2{k))dk. Here, Xn is the magnetic suscep- 
tibility in the normal state, and Y{k,T) is the angle de- 
pendent Yosida function, 

n,,T) = iyi-^)\^^^^^, (32) 

where f{E) is the Fermi distribution function. It should 
be noticed that the magnetic susceptibility decreases for 
the magnetic field parallel to the rf-vcctor. For the fxy+- 
state, it is reasonable to assume that d = p^x + PyV or 
d = PyX — PxV is stabilized below Tc among any linear 
combinations of these states. This is because the conden- 
sation energy is maximally gained in these states within 
the weak coupling theory. Then, the susceptibility along 
the plane Xab decreases to the half of its value in the nor- 
mal state, as shown in Fig. 7. In order to obtain Fig. 7, 
we have assumed the order parameter below Tc as d{k) = 
A{T){(f)^(k)x + (j)y{k)y) where (t)^{k) and (f>y{k) arc ob- 
tained in eq. (21) as A^i-(A;, IttTc) = —(j)x{k) + i(/)y(fc). 
The wave functions ^x(fe) and (j>y{k) have the symmetry 
of Px-wavc and py-wave, respectively. We have calciilatcd 
the temperature dependence of A(T) by using the effec- 
tive model having the separable pairing interaction and 
solving it within the BCS theory. The same procedure 
has been used in Ref. 59. For the i^xy-state, the in-plane 
magnetic field favors one of the degenerate states. For 
example, d = af2X + fiy is favored by the magnetic field 
along X-axis. Then, Xah decreases owing to the induced 
/2-wave component, however the decrease is very small 
as shown in Fig. 7 because \a\ <^ 1. Here, we have cal- 
culated the temperature dependence of order parameter 
in the same way as for the Pxy-i--state. The susceptibil- 
ity along c-axis Xc does not decrease in both -Pxy+- and 
Fxy-states. Although Xc decreases to zero in the F^-state 
at T = 0, the observation of this decrease will be diffi- 
cult because the (i- vector rotates in a very weak magnetic 
field. 

The Knight shift below Tc has been measured by some 
groups.*''^' Most of the measurements have been 
performed under the magnetic field parallel to the plane. 
It has been reported that the Knight shift in Co-NMR 
and 0-NMR decreases below Tc in a weak magnetic 
field.^'^'^'^'^"'^^ This observation is consistent with the 
Pxy-i--state in our results. The Pxy-i--state is furthermore 
consistent with the /iSR measurement"'^''^ which does not 
detect any time-reversal symmetry breaking expected in 
the Pz-state and d^2_y2 + irfxy-state. Contrary to the 
other NMR measurements, the Knight shift in D-NMR 
is almost temperature independent below Tc.^-* Although 
this observation is incompatible with the other NMR 
measurements, this is consistent with Fxy- and F^-states 




T/T 

c 

Fig. 7. Temperature dependence of magnetic susceptibility along 
the plane Xab divided by the normal state value Xn. We consider 

the magnetic field along x-axis and assume d = p^x + Pyij in the 
Pxy+-state and d = af2X + fiy in the Fxy-state. We obtain the 
Fermi surface and momentum dependence of order parameter for 
the parameters A and B in Fig. 4, respectively. 



in our calculation. 

4-4 Crossover from X <^W to X:$> W 

In §4.1-3. we have considered the spin-orbit coupling 
much smaller than the band width, namely A <C W . This 
relation is expected in most of the d-electron systems in- 
cluding NaxCo02 • 2;H20. On the other hand, the situa- 
tion is opposite in the /-electron systems which include 
odd-parity superconductors. It is not clear whether the 
role of spin-orbit coupling is qualitatively different be- 
tween the weak coupling region X <ti W and strong cou- 
pling region X ^ W. Therefore, it is interesting to in- 
vestigate global behaviors of two-orbital Hubbard model 
from X <^ W to X ^ W, although the strong coupling 
region X^ W is unreahstic for NaxCo02 • yH20. 

In Fig. 8, we show the anisotropy of eigenvalues of 
linearized Dyson-Gorkov equation Aan = The 
small value of Aan means that the splitting of supercon- 
ducting Tc is large. It is shown that the splitting of Tc 
takes its maximum for an intermediate value of A. Thus, 
the anisotropy of d-vector decreases with increasing the 
spin-orbit coupling in the strong coupling region X^ W 
in contrast to the weak coupling region. This behavior is 
quite different from the anisotropy of spin susceptibility 
Xan = Xc/Xab which has been shown in Fig. 8 for a com- 
parison. The anisotropy of spin susceptibility is enhanced 
monotonically with increasing the spin-orbit coupling, as 
expected. Note that the structure around A = 2 is due 
to the rapid change of Fermi surface. There are two hole 
pocket Fermi surfaces enclosing the K-point for A > 2. 

It is easy to understand the behavior of Aan by taking 
the limit A — > 00 in advance. In this limit, the two-orbital 
Hubbard model is reduced to the single-orbital Hubbard 
model written as, 

Hi = E^2(fe)cJ.^^Cfc_^ + ^^±^ ^n^.tna, (33) 
k,s ' ' 

where the Coulomb interaction is renormalized to '^^^ . 
It is notable that s in eq. (33) is not the spin but the 
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Fig. 8. The anisotropy of eigenvalues of linearized Dyson-Gorkov 
equation, Aan = aJ^/aJ^ in a wide range of spin-orbit coupling. 
The solid (dashed) line shows the result for the p-wave (/-wave) 
symmetry. The parameters are Jyi/U = 7/24, ec = 0, {Ju/U = 
5/24, ec = -0.5) and T = 0.01. Here, we fix the crystal field 
splitting Be and total hole density n = 0.5 instead of Ue- The 
anisotropy of spin susceptibility Xan = Xab/Xc is shown by the 
thin solid line. The double dotted line shows A|/aI^ where A| is 
the eigenvalue in the single-orbital Hubbard model (eq. (33)). 



pseudo-spin. This single-orbital model has the SU{2) 

symmetry for the rotation of pseudo-spin. Therefore, spin 
triplet superconducting state is isotropic for the direc- 
tion of drvectoT defined by the pseudo-spin. In Fig. 8, we 
see that Aan in the two-orbital model comes to unity in 
the limit A — > oo. Then, the maximum eigenvalue of lin- 
earized Dyson-Gorkov equation actually comes to that 
of single-orbital Hubbard model. On the other hand, the 
^-factor of Zeeman coupling term is anisotropic for the 
pseudo-spin, and therefore the anisotropy of spin sus- 
ceptibility remains in the single-orbital model. In other 
words, the effect of spin-orbit coupling is absorbed by the 
character of quasi-particles and the residual interaction 
is isotropic in the strong coupling limit A ^ oo. 

Although the two-orbital Hubbard model in eq. (11) is 
too simple to provide a general formula, the same situa- 
tion is generally expected when the heavy quasi-particles 
consist of one Kramers doublet. Actually, it has been re- 
ported by the NMR measurements^^) that the d-vector 
in UPt3 is almost isotropic, although the spin suscepti- 
bility is significantly anisotropic. Thus, the anisotropy of 
d-vector and that of spin susceptibility is quite different 
in general. 

5. Comparison with Sr2Ru04 

Interestingly, we have found some similar aspects be- 
tween NaxCo02 • yH20 and Sr2Ru04. The latter is 
the most established spin triplet superconductor in d- 
electron systems. ''^ One of the similarity is the orbital 
dependent superconductivity. The three t2g-orbitals is 
divided to two groups, namely dxy-orbital and dy^- 
dzx-orbitals in Sr2Ru04, oig-orbital and Cg-orbitals in 
NaxCo02 • J/H2O. The superconductivity is mainly in- 
duced by the former in Sr2Ru04 while by the latter in 
NaxCo02 • yH20. Another interesting finding is that the 



spin-orbit coupling term has the same matrix element as 
eq. (10). This enables us to summarize the results on the 
d-vector in a unified way as shown in Table III. 



Sr2Ru04 


NaxCo02 •2/H2O 


Square lattice 


Triangular lattice 


(ixy 




< 






iz 


p-wave 


p-wave 


/-wave 


O(A^) 


0(A) 


0(A) 


O(A-) 


d ||z 


d l|xy 


d l|xy 


both 



Table lU. An unified description on the d-vector between 
Sr2Ru04 and Na^xCo02 • J/H2O. The second column shows the 
symmetry of crystal. The third column shows the orbitals leading 
to the superconductivity. The forth column shows the important 

component of orbital moment. The fifth column shows the sym- 
metry of superconductivity. The sixth column shows the leading 
order term with respect to the spin-orbit coupling. The last col- 
umn shows the direction of d-vector. The results on Sr2Ru04 
have been obtained in Ref. 50. 



We also show the case where the dy^- and d^x-orbitals 
induce the superconductivity in Sr2Ru04. In this case, 
the results on the d-vector are qualitatively the same as 
the p-wave superconductivity in NaxCo02 • yH20. Actu- 
ally, the analysis of the first order term in A (see §4.2) can 
be applied to the case in Sr2Ru04 in the same way.^°) In 
both cases, the orbital moment along the c-axis stabilizes 
the d-vector along the plane. 

The first order term in A vanishes in Sr2Ru04 when 
the dxy-orbital is active. Although this is similar to the 
/-wave state in NaxCo02 • yH20, the microscopic ori- 
gin is quite different. The first order term disappears in 
Sr2Ru04 since the hybridization term with dy^- and d^x- 
orbitals vanishes. 

Note that the results on Sr2Ru04 are consistent with 
experiments. It is believed that the dxy-orbital is mainly 
superconducting in Sr2Ru04. Then, we obtain the chiral 
superconducting state, ^'^^ namely Pz-state in Table II, 
which is consistent with the /iSR measurement.*'^) Owing 
to the disappearance of the first order term in A, the d- 
vector can be rotated by a weak magnetic field along the 
c-axis. Actually, the NMR measurement has detected a 
d-vector along the plane under the weak magnetic field 
H II 

6. Summary and Discussions 

In this paper, we have investigated the d-vector in the 
possible spin triplet superconductor NaxCo02 • J/H2O on 
the basis of the two-orbital Hubbard model representing 
the eg-orbitals. There remains a 6-fold (3-fold) degener- 
acy in the p-wave (/-wave) superconducting state if we 
neglect the spin-orbit coupling. Therefore, we include the 
L-S coupling term in Co ions into the Hamiltonian and 
determine the pairing state on the basis of the second 
order perturbation theory. 

We find that the role of spin-orbit coupling is quite 
different between the p-wave superconductivity and /- 
wave superconductivity. The d-vector is always along the 
plane if the orbital part has the p-wave symmetry. On 
the other hand, the direction of d-vector depends on the 
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parameters if the orbital part has the /-wave symmetry. 
In our calculation, the superconducting Tc is determined 
for each direction of rf-vector. We show that the splitting 
of Tc in the j;-wave state is much larger than that in the 
/-wave state. 

Such a different role of spin-orbit coupling is explained 
by analyzing the first order term with respect to the spin- 
orbit coupling A. This term is effective in the p-wave 
symmetry, but ineffective in the /-wave symmetry. This 
property is hold in all of the perturbation terms with 
respect to the Coulomb interactions. Therefore, the re- 
sults in this paper will be valid beyond the second order 
perturbation theory adopted here. The only assumption 
is that the superconductivity is mainly induced by the 
Cg-orbitals. 

The determination of d-vector is especially important 
for the interpretation of Knight shift measurements be- 
cause the magnetic properties of spin triplet supercon- 
ductor depend on the direction of rf-vector. Our results 
indicate that the d- vector can be rotated by a weak mag- 
netic field in the case of /-wave superconductivity, while 
the rf-vector is strongly fixed against the applied field in 
the case of j)-wave superconductivity. The NMR Knight 
shift along the plane will decrease in the latter case be- 
caiise the rf-vector has both x and y components. On the 
other hand, the Knight shift will be almost temperature- 
independent in the former case in all directions of ap- 
plied magnetic field. Unfortunately, experimental results 
seem to be confusing. The NMR Knight shift in Co-site 
6, 7, 60, 61) g^jjjj 0-site has shown a sizable decrease for 
the parallel magnetic field. This is consistent with our re- 
sults for the p-wave state. On the other hand, the Knight 
shift in D-site has reported a qualitatively different re- 
sult^^ which is consistent with our results for the /-wave 
state. 

Finally, we suggest some future experimental studies 
which are highly desired. First, the NMR Knight shift 
along c-axis may provide conclusive evidence for the pair- 
ing symmetry. Kobayashi et al. has reported a pioneering 
result which shows a slight decrease of Knight shift be- 
low TcJ^ Although this observation implies the spin sin- 
glet superconducting state, the theoretical interpretation 
seems not to be conclusive. It is noted that the decrease 
of c-axis Knight shift in Ref. 7 is very small compared 
with that of in-plane Knight shift. Although this small 
decrease may be explained on the basis of the vortex state 
with spin singlet pairing, another theoretical interpre- 
tation may be also valid. For example, we suggest three 
possibilities based on the spin triplet superconducting 
state. The first one is the coexistence of p- and /-wave 
superconductivities which will be discussed in the subse- 
quent paper. ■"'^^ Extending the present theory to deter- 
mine the pairing state below we find this coexistent 
state around the boundary between the p- and /-wave 
states in Fig. 4. The magnetic susceptibility decreases in 
this state for all directions of magnetic field. The second 
one is the role of strong electron correlation. The con- 
clusions obtained in the weak coupling BCS theory can 
be modified in the strongly correlated electron systems, 
especially near the magnetic instability. The third one 
is the role of disorder which leads to the rotation of d- 
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vector. The detailed measurements of c-axis Knight shift 
including the magnetic field dependence, sample depen- 
dence and impurity effect are highly desired for a clear 
identification of the pairing state. 

Second, the search for the multiple phase transition 
under the magnetic field is particularly interesting in the 
future study. The phase transition accompanied by a ro- 
tation of d-vector is expected in the spin triplet super- 
conductor under the magnetic field. Such a phase tran- 
sition is promising when the iJc2 if' large, namely when 
the magnetic field is applied along the plane. We have 
actually determined the phase diagram under the par- 
allel magnetic field on the basis of the weak coupling 
theory. We find that such a phase transition is ex- 
pected in the p-wave superconducting state and p + f 
coexistent state while that is not expected in the /-wave 
state. Interestingly, the second phase in the j>wave state 
is different from all of the pairing states in Table II, and 
the NMR Knight shift decreases even in this phase. The 
results on this subject will be reported in another publi- 
cation.^^' 

Third, the precise measurement on the; critical field 
Hc2 along the plane is important. Although some authors 
have reported the measurement, the results are contro- 
versial.^^' If the Hc2 far exceeds the Pauli paramag- 
netic limit, the spin triplet pairing state is promising. 
This issue is closely related to the second issue, namely 
the possibility of multiple phase transition. For exam- 
ple, it has been reported that the Knight shift in the 
Co-NMR and fiSR does not decrease in a high magnetic 
field along the plane. If this magnetic field is still 
below Hc2, the qualitatively different behavior from the 
low magnetic field region indicates a multiple phase tran- 
sition. 

At last, it is desirable to determine the Fermi surface of 
superconducting material experimentally. Although the 
existence of Cg-Fermi smfaee is a basic assumption of this 
paper, it is not clear whether the Cg-Fermi surface exists 
or not. Recent ARPES measurements^^' ^'^^ on the non- 
superconducting material have reported a qualitatively 
different Fermi surface from the LDA calculations. ^'*"^^' 
We think that the measurement on the superconducting 
material is highly desired. 
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